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SYZYGIES IN THE TWO CENTER PROBLEM 


HOLGER R. DULLIN AND RICHARD MONTGOMERY 


Abstract. We give a complete symbolic dynamics description of the dynamics 
of Euler’s problem of two fixed centers. By analogy with the 3-body problem we 
use the collinearities (or syzygies) of the three bodies as symbols. We show that 
motion without collision on regular tori of the regularised integrable system are 
given by so called Sturmian sequences. Sturmian sequences were introduced by 
Morse and Hedlund in 1940. Our main theorem is that the periodic Sturmian 
sequences are in one to one correspondence with the periodic orbits of the two 
center problem. Similarly, finite Sturmian sequences correspond to collision- 
collision orbits. 


1. Introduction 

Associating a symbolic dynamics to a smooth flow goes back at least to Hadamard 
[8] with serious contributions due to Morse and Hedlund [13, 14]. The advent of 
Smale’s Horseshoe [16] and Anosov flows [1] made symbolic dynamics into an area 
of study connected to smooth dynamical systems in a central way. In the area 
of celestial mechanics, a number of partial results have been obtained about a 
tentative symbolic dynamics for the planar three-body problem. See for example 
[11], [12], and [10]. The symbols of this symbolic dynamics are the “letters” 1, 
2 and 3 with letter i representing an instant during the motion where the three 
bodies have become collinear with body i in the middle. For historical reasons we 
call the resulting sequences “syzygy sequences”. A big open question is, for the 
planar Newtonian three-body problem “what are the possible syzygy sequences 
realized by motions of the three bodies?” Here, we answer this question for a 
much simpler integrable system: the planar two-center problem as investigated by 
Euler [7]. Classical reference on this system are Charlier [4] and Whittaker [18], 
and a modern treatment is given in [17]. The symbols here have precisely the same 
meaning as for the full planar three-body problem but the system is completely 
integrable, hence we can completely solve the question. We hope that this full 
solution in the integrable limit may shed light on the honest three-body problem. 

2. Symbolic dynamics 

Consider a flow M —)■ M on a manifold M. Fix a finite collection of co¬ 
dimension one surfaces Si C M, i = 1, ... ,?7, (possibly with boundary) which we 
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call windows. Attach symbol Sj to window Si. Associate a finite symbol sequence 
to each finite orbit segment by recording the windows in the order in which they 
are hit. Thus orbit segment 0*(a:o), 0 < t < T has symbol sequence Sjj,... 
provided 0*“(a;o) € a = 1,..., A^ where Q <ti <t 2 < .. .tjq <T are the times 
at which the orbit hits some window. (The sequence is empty if no window is hit.) 
This finite symbol sequence is unique provided the orbit avoids the intersections 
of the windows. The sequence is stable with respect to small deformations of the 
initial condition xq provided the orbit (p^ixo) is transverse to the windows. Letting 
T —>■ oo generally gives an infinite sequence. Running time backwards gives a 
bi-infinite sequence of symbols. 




Figure 1. Left: Linear flow on the torus with initial condition 
(0,0.15) and slope W = l/vr. The symbol sequence of the shown 
orbit segment is VVVHVV. Right: Corresponding Poincare map 
from a: = 0tox = l = 0 giving the exponents n* = 0, 0,1, 0. 


2.1. Symbolic dynamics for integrable systems. Suppose now that our flow 
on M is a Liouville integrable system with 2 degrees of freedom. Then M is a 
symplectic manifold M of dimension 4 endowed with 2 smooth functions G, H : 
M —)■ M in involution, independent almost everywhere, one of which, H, we take 
to be the energy of the system. The pair F = {G, Lf) : M —>■ is called the 

integral map. The Liouville-Arnold theorem states that if c is a regular value 
of the integral map, then every compact connected component of the pre image 
F~^{c) is a 2-torus. Denote one of these tori by Tc. Then in a neighbourhood U of 
Tc one can introduce action-angle variables [Ii, 12 , 61 , 62 ), such that the dynamics 
is given by 


~ -^ 2 ), Ij — 0 , 


J = l,2. 
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The tori are labelled by the value of the action I = (/i,/ 2 ). On each torus the 
dynamics is straight line motion when viewed from the covering space of that 
torus: 6 i{t) = 6i{0) +Ui{I)t. When changing from one torus to another, I and in 
general uji{l) will change, and hence the slope of the lines changes. 

If the windows Si G M intersect the torus transversally then they form a 
system of windows Tc 0 Si (possibly empty) for the restricted linear dynamics 
on Tc- In our situation of the planar two-center problem the windows will be 
horizontal or vertical lines: 6 i = const or 62 = const. With this in mind, the 
following example which goes back to Morse and Hedlund [14] is central to our 
work. 

Example: Consider a linear flow on mod with positive slope 

m = uj 2 /u:i. As windows take the two two circles obtained by projecting the x-axis 
and the y-axis of onto the torus. Denote their symbols ‘H’ for horizontal and 
‘V’ for vertical. Lifted to the xy plane, the windows yield the lines of graph paper 
intersecting at the lattice points. The lifted orbits are the lines 

(1) y = mx -1- h. 

Every such orbit has a unique symbol sequence as long as it stays away from lattice 
points. 

The ratio of Lf’s to E’s in the symbol sequence of an orbit tends to the slope 
m as the length of the orbit tends to inhnity. The precise sequence of H's and 
Id’s depends on the initial condition t/q, i.e., on b. The resulting sequences are 
called “Sturnian sequences” after Morse and Hedlund coined this term in [14]. To 
describe this sequence let [y\ G Z denote the greatest integer less than or equal 
to ?/, for y a real number. A bit of experimentation drawing lines on graph paper 
suggests that the associated symbol sequence always has the form 


( 2 ) ...VH^^VH^WH^K.. 

where each nonnegative integer rij is either [mj or [mj -|- 1. Which one though? 
To describe the n,, let {y} denote the fractional part of a real number y so that 
y = \_y\ +{y} with 0 < {y} < 1. Use the Poincare map associated to cutting across 
the vertical circle: y ^ y + m so as to obtain a sequence of points yi mod 1 on 
the circle, these being the orbit of j/o uder the Poincare map, and their associated 
fractional parts, {yk} = {yo -l- km}. Then 


( 3 ) 


nk 


[mJ, if {yk-i} + {m} < 1 
[m}+l, if {m} > 1 


Here we have started with the point at (0, yo) and marked the initial ‘U’ in equation 
(2) as corresponding to “sq” accordingly. An equivalent formula for the exponents 
is nk= [m + 


Definition 1. The sequence of integers {ufcjfcez defined by equation (3), or the 
corresponding symbol sequence of H’s and V’s will be called the “Sturmian sequence 
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Figure 2. Orbits in the (x, i/)-plane for the symmetric case mi = 
m 2 = 1/2, (i = 1, and h = —0.23. Type Satellite (left) intersects 
the T-axis near one of the centers only. Type Lemniscate (middle) 
intersects the a:-axis near both centers. Type Planetary (right) does 
not intersect the segment between the centers. 


for rotation number m and intercept yo ”, assuming that the line y = mx + i/o does 
not hit any lattice points. 

Lemma 1. Suppose that the slope m is rational. Then when viewed as a periodic 
word, the Sturmian sequence associated to rotation number m and intercept yo is 
independent of the choice of intercept. Thus we can speak of “the Sturmian word 
for m” in case m is rational. 

Proof. Suppose that m = p/q is rational, with p,q relatively prime. Then 
the initial conditions (0, 6), 0 < 6 < 1 which pass thru the lattice points are those 
for which b = i/q,i = 1,2 ,... ,q — 1 and they cut the basic interval into q equal 
intervals. Within any one of these intervals the sequence is constant, since the 
associated line y = mx + b can be translated without hitting lattice points. What 
about when we move our initial condition (0, b) to a different interval? Use the 
lattice translations to decompose all the vertical unit intervals between vertically 
adjacent lattice points into these q equal segments. Because p and q are relatively 
prime, any line starting in the interval 0 < y < 1 /q eventually hits all of the q 
other segments. We can thus imagine starting at any one of these q segments, but 
at a later time along the line, so arriving at a shifted version of the same sequence 
which corresponds to b lying in this other segment. QED 

Remark. When m is irrational the Sturmian word depends on the intercept y^ in 
an interesting and nontrivial manner. 

Remark. Sturmian sequences arise in billiards [5], in the Henon map [6], and in 
twist maps of the cylinder [2]. 
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3. The Two Center Problem and its tori. 


We consider a unit test particle moving in the xy plane attracted by two 
fixed masses of mass mi and m 2 called “centers” and located at the points {—d, 0) 
and (d,0). The Hamiltonian in canonical coordinates x,y,Px,Py with symplectic 
form dx A dp^ + dy /\ dpy is 

/ A\ TT ^ / 2 I 2\ Tfli Tfl2 

(4) H=-{p^+p)- ^-_ ===- — 

^ V(a; + d)2 + j/2 ^(^x-dy + y^ 

Euler [7] proved that the system is integrable with independent second integral 


(5) G =]^{xpy-ypyy+ ]^d^pl + dx 


mi 


m2 


\/{x + dy + y"^ \/{x — dy + y‘^ 


{G, H) form our “integral map”, a map from phase space to the space of values of 
these two integrals. 

Integration involves switching to regularizing variables for configuration space 
which are denoted A, u below. For details see section 7 below. These variables 
have the local effect of a Levi-Civita regularization about each of the two centers. 


Definition 2. By a regular torus T we mean a compact connected component of 
an inverse image of a regular value of the integral map, expressed in regularizing 
variables. By a regular periodic orbit we mean a periodic orbit without collisions 
lying on a regular torus. 


Excluding the critical values of the integral map, there are precisely three topo¬ 
logical types of tori, denoted by S, L and P. The projection of a torus of each 
type onto configuration space are as shown in figure 2. After regularization we can 
define “action-angle” variables with angles di, 62 . These angles identify a regular 
torus with We will say a torus, expressed in these variables, has been 

“flattened”. 


3.1. Windows, Syzygies. The phase space is T*(]R^\ A) where A = {(—d, 0), (d, 0)} 
are the collisions. When the test particle crosses the x-axis the three masses be¬ 
come collinear. The centers divide the x-axis into three intervals, labelled “1” for 
the right infinite interval (—C) 0 , d), “2” for for the left infinite interval (d, 00 ). and 
“3” for the central finite interval —d < x < d. This labelling corresponds precisely 
to the syzygy sequence labelling of the references [12] and [10] mentioned above 
provided the test particle is labelled “3”, the right center as “1” and left center as 
“ 2 ”. 

The inverse image of the collinear locus, with its three windows, drawn on the 
flattened torus is shown for each type in figure 3. These pictures are at the heart 
of our analysis. We continue to denote the inverse image of our syzygy windows 1, 

2, and 3, under the regularization map, intersected with the torus, and flattened 
by the symbols 1, 2, or 3. These windows are a collection of orthogonal lines. 
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Figure 3. The toms (dashed) in angle variables with a tra¬ 
jectory and windows with corresponding symbols for type L (Lem- 
niscate), S, (Satellite), P (Planetary). The corresponding symbol 
sequences for the orbit segment shown (line with an arrow, contin¬ 
ued with periodic boundary conditions) are 312132, 311131, 1212, 
respectively. 


How did a decomposition of a straight line - the x-axis - into intervals become 
a system of orthogonal lines? A regularization, qualitatively similar to the the 
Levi-Civita regularization, is needed to write down the action angle variables and 
hence express the torus. At the core of Levi-Civita regularization is the squaring 
map, the conformal map ii] ^ {x ^ iy) = (.^ -f- ir])"^ under which the inverse 
image of a straight line through the origin becomes a “cross”: a pair of orthogonal 
lines, and this fact explains how the x-axis turned into a collection of orthogonal 
lines. See section 7 for details. 

4. Main results: syzygy sequences for the two center problem 
The main result is 

Theorem 1. The periodic syzygy sequences for regular periodic orbits of the two 
center problem are of the form: 

L type: 13”i23"'M3”3 2 ... 3""* (Ts and 2’s alternating) 

S type: or 23^^23^^23^^.. .23^‘ 

P type: (12)^ 

where the exponents Ui are those of the Sturmian sequence (3) associated to the 
rational rotation number W for that orbit’s torus (Definition 1 ). The length of the 
fundamental sequence is 2{p + q) where W = p/g > 0. (For the range of possible 
rotation numbers see the final sentence of the next theorem.) 
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W=1/5 

:213212123121; 


W=1/3 

:23121321; 






W=3/1 

;33323331: 


W=5/1 

:333332333331: 


W=2/1 

;332331: 



W=4/1 

;3333233331; 


W=6/1 

;33333323333331; 







W=5/2 

;31332333133233; 



Figure 4. Periodic orbits of Lemniscate type in the {x, y) plane for 
the equal mass case mi = m 2 = 1/2, d = 1, h = —0.23 with various 
rotation numbers W. To the right more oscillations along the |/-axis 
are added, to the left more oscillations along the a:-axis are added, 
keeping the parity of W the same along each row. 


See hgure 4 for examples of periodic type L orbits, their rotation numbers and 
syzygy sequences for equal masses. See hgure 5 for examples of periodic type S 
orbits, their rotation numbers and syzygy sequences in a case of distinct masses. 

Missing periodig orbits. The only periodic orbits whose sequences are not 
represented in the theorem above are those on non-regular level sets of the integral 
map. There are hve such families of periodic orbits corresponding to critical values 
of the integral map, see hgure 11. Three of these hve correspond to periodic orbits 
that are collinear, i.e. they move along the axis y = 0 that contains the two centers, 
and as such don’t have well dehne syzygy sequences. The other two families have 
either syzygy sequence 3 (separating type L from type S, orbit is on the symmetry 
axis X = 0 when then masses are equal), or syzygy sequence 12 repeated (orbit 
encircling both centers, the outer boundary of type P). 

Our main result is a corollary of the next theorem which includes the inhnite 
aperiodic orbits - the dense windings on the tori. 

Theorem 2. The syzygy sequences realized by regular non-collision orbits for the 
two-center problem are precisely one of the following possibilities: 

... 121212 ... ,in the P case, 
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Figure 5. Periodic orbits of Satellite type in the {x,y) plane for 
the asymmetric case mi = 1/3, m 2 = 2/3, d = 1, h = —0.23 with 
various rotation numbers W. Parity if the same along each row. The 
third row shows tori of type SI, all other of type S2. 

... 13^^23^n3^^2 in the L case, 

... ..., or ... 23"i23”223"32 ... in the S case 

In the last two cases the integer sequence Ui occurring is the Sturmian sequence 
(Definition 1) associated to the rotation number W = W{g, h) of that orbit’s torus 
T{g,h), and any possible intercepts yo. The range of rotation numbers W deter¬ 
mining the Ui can be read off from tables 10 and 11 below. 

Remark. The rotation number W only depends on the regular values g, h and 
not on the choice of torus when F~^{g, h) consists of more than one torus. This 
“coincidence” is a consequence of the fact that the two-center problem is integrable 
in elliptic functions, and moreover that the differentials that dehne the rotation 
numbers are of hrst kind. 

The main theorem has some simple consequences. 
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ea 



Figure 6. Two solution curves with W = 2/?) bounding a region 
of orbits without collisions. Half-lattice points correspond to colli¬ 
sions. 

Corollary 3. In any orbit on a regular torus of type L or S the number of con¬ 
secutive 3s is either n or n + 1 where n is the integer part of W. 

Proof. This follows directly from the dehnition of the Sturmian sequence. □ 

Corollary 4. In all cases symbol 1 and 2 never appear adjacent to each other (“no 
1 or 2 stutters”). Lemniscate orbits with W < 1 have no stutters, while for W > 1 
they have 3-stutters. Every Satellite orbit has at least one 3-stutter. 

Proof. For type L the symbols 1 and 2 alternate, so there cannot be stutters of 
these symbols. For type S there would be 1- or 2-stutters if any Uj = 0, however 
we have numerically verihed that W > 1 for type S, so by equation (3) we have 
Hi > 0. At least one 3-stutter, i.e. one Uj = 2, must occur when W > 1, because 
then the line in hgure 3 with slope > 1 must intersect the vertical lines of symbol 
3 at least once. □ 


5. Guts of the proof. 

Each regularized Liouville torus T = T( 5 f, h) maps to conhguration space by 
composing its inclusion into regularized A, u, p\, pu-ph.a.se space and projecting the 
regularized phase space onto the usual xy conhguration space. The windows on 
T are the inverse image of the collinear locus (the x-axis) with respect to this 
map. We will show in sections 7 and 8 that these windows are a hnite collection 
of horizontal and vertical lines, positioned as indicated in hgure 3 according to the 
cases S, L or P. Here “horizontal” or “vertical” mean that when we coordinatize 
T in standard angle variables 6162 of “action-angle” as then these lines are 

parallel to the 61 or the 62 axes. Lifted to the universal cover, then we get almost 
exactly the ‘graph paper picture’ used for writing Sturmian words. In the S and L 
cases we have 4 lines in all, a horizontal pair and a vertical pair, see hgure 3. The 
parallel lines in each pair are separated by 1/2 a lattice unit. In the P case we 
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Figure 7. Orbits of type L with W = 1 and syzygy sequence 1323 
for h = —0.23 and mi = m 2 = 1/2 for a sequence of initial conditions 
with different vq = zvr/ld, i = 0,..., 7 and the same Aq = 0 on the 
same torus. For z/q = 0 the orbit has additional discrete symmetry, 
while for z/q = vr/2 it is a collision-collision orbit. 


have only the two vertical lines separated by 1/2 a unit. The intersections of the 
lines represent collisions of our test particle with the appropriate primary. This 
pattern is represented on the plane as a doubly periodic pattern of parallel lines. 

On a fixed torus T = T( 5 f, h) the solution curves are a family of parallel straight 
lines y = Wx + b with W = W (T) constant, and b varying. Here we have changed 
variables so x = 9i, y = 62 so as to conform to the discusion of subsection 2.1. 
The torus supports periodic solutions if and only if W is rational. A hnite number 
of these solutions will have collisions. Those remaining have a syzygy sequence. 

The following is essentially a restatement of theorem 2: 

Lemma 2. On a fixed regular rational torus T = T( 5 f, h) all eollision-free solutions 
have the same periodic syzygy seguence. This seguence is as described in theorem 2, 
where W = W{g,h) is the rotation number ofT. If W = p/q then the length of 
this common syzygy sequence is 2{p + q) in the S and L cases, and 2q in the P 
case. 

Proof. Let W = p/q he given, and the labelled set of vertical and (possibly) 
horizontal lines be given. Each solution is a straight line of slope p/q, see figure 6 

P case. There are no horizontal lines, so no collisions. Vertical lines marked 
with 1 and 2 alternate. The only possible sequence is (12)® for some power s. Each 
straight line cuts across q basic units in the horizontal direction before closing, so 
that we have s = q. 

S and L cases. Fix attention on one of the 4 collision points: these being the 
intersection of the vertical and horizontal lines. Shift the fundamental domain 
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describing the torus so that this point is at the origin. Collision solutions are 
those passing through a lattice point and we throw these out. Because the line 
pairs are half a unit a part, it makes sense to magnify the lattice. Consider the new 
lattice ((1/2)Z)^ in the plane, the torus’s universal cover, whose points represent 
collisions. We are now in a situation identical to that of subsection 2.1, except 
that the lines have different labellings! The independence of the sequence on the 
initial condition of the torus follows from lemma 1 of subsection 2.1. Keeping 
track of the labelling leads to the sequences described in theorem 2. Everything 
has already been done by Morse and Hedlund, as described in section 2.1, with the 
exception of the counting of sequence length. For this count, observe that there are 
2 horizontal lines per vertical unit and 2 vertical lines per horizontal unit. Every 
solution must traverse q horizontal units and p vertical units and so crosses 2q + 2p 
lines in all. 

□ 

Remark. The fact that the symbol sequences are the same for solutions on a fixed 
rational torus does not mean that these solutions are the same! Indeed, they sweep 
out a torus-full of solutions. In figure 7 we show the deformation of the type L 
orbit with W = 1/1 when the initial condition on the torus is changed. The orbit 
undergoes a deformation from a symmetric state all the way to a collision orbit. 
Note that for a generic non-integrable perturbation one expects that almost all of 
these periodic orbits will be destroyed. 

6. Collision orbits 

In figure 8 we show some collision orbits, which are found on tori with rational 
W for initial conditions with A = 0 and v = ±7r/2, as we now show. 

Corollary 5. Orbits that connect two collisions lie on a regular torus with rational 
rotation number W. 

Proof. Collision points are half integer lattice points in the covering plane of the 
torus. Any line connecting such lattice points has rational slope. □ 

Corollary 6. Orbits that have a single collision lie on a regular torus with irra¬ 
tional rotation number W. 

Proof. A line with rational slope that hits a half integer lattice collision point 
in forward time, will also hit another such point in backward time. Hence only 
lines with irrational slope can hit one lattice point (in either forward or backward 
time). □ 

Corollary 7. Regular orbits that keep a finite distance from collisions lie on a 
regular torus with rational rotation number W (unless in the P family, which has 
no collisions). 
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W=7/3 



Figure 8. Collision orbits of Lemniscate type in the {x, y) plane for 
the equal mass case mi = m 2 = 1/2, d = 1, h = —0.23 with various 
rotation numbers W. To the right more oscillations along the y-axis 
are added, to the left more oscillations along the x-axis are added, 
keeping the parity of W the same along each row. Collision orbits 
with W = odd/even visit the same collision twice. The last row 
shows some brake-collision orbits of Satellite type (same parame¬ 
ters). 

Proof. A line with irrational slope will come arbitrarily close to a half integer 
lattice point. Hence the only lines that keep a hnite distance have rational slope, 
unless they are collision-collision orbits, as described in the previous corollary. □ 

Corollary 8. Any regular torus of type S or type L with rational rotation number 
has exactly two collision-collision orbits. 

Proof. This follows from lemma 2, the fundamental region covered by two copies of 
the strip shown in hgure 6 has exactly two collision-collision orbits, corresponding 
to the two lines shown in the hgure. □ 
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Figure 9. Level lines of Hx (left) and Hy (right) with critical points 
marked by red and blue dots, respectively, and pre-images of critical 
values in the corresponding color. Green contours denote the levels 
corresponding to the three example shown in Fig. 1. Always mi = 
m 2 = 1/2, d = 1. Top: h = —1/6 > h\, types S, L, and P occur. 
Middle: h = —2/3 > h*, types S, and L occur. Bottom: h = —6/5 < 
h*, type S only. 


7. Euler’s Problem of Two Fixed Centers 

In this section we recall what we need from the detailed description of the 
solution to the problem of two hxed centers given by [17]. Many of the results here 
are refered to in theorem 2. 

Without loss of generality we set d = 1 so that the centres are place at (1,0) 
and (-1,0). Standard confocal elliptic coordinates (A, z/) are then dehned by the 
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-2-1 0 1 2 -3-2-10123 



Figure 10. Level lines of Hx (left) and (right) with critical 
points marked by red and blue dots, respectively, and pre-images 
of critical values in the corresponding color. Always mi = 3/10, 
m 2 = 7/10, d = 1. Top: h = —1/6 > h^, types S’, L, and P 
occur. Middle: h = —2/3 > h*, types S', S, and L occur. Bottom: 
h = —6/5 < h*, type S and S’ only. 


relation x + iy = cosh(A + iz/). We will slightly alter the standard relation, dehning 
instead 

X + iy = sin(z/ -|- iX). 

In terms of real and imaginary parts we have 


{x, y) = (cosh A sin z/, sinh A cos z/) 
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Combined with the time reparameterization r = r(t; A, u) given by dt = (cosh^(A) — 
cos^(z/))dr, these variables regularize and separate the Hamiltonian at energy h: 

H = H{X,iy,px,Pu;h) := {H - h)^ = Hx{X,Px) + H^{u,p^) 

dr 

into two one-degree of freedom Hamiltonians: 

Hx{p, q',h) = —— (rui -|- m 2 ) coshg — h cosh^ g, 

H^{p, q; h) = — + (mi — m 2 ) sing -|- hsin^ g . 

We must take the value of if to be 0 since we want H = h. Then we have that 
H\{p,q) = —g and Hu{p,q) = g where g is the value of the separation constant 
which is also the second integral G, see equation (5). (Our Hx,H^ differ from 
those in [17] by —h.cosh^(A), hsin^(p) respectively.) 

The map (A, z/) —)■ {x,y) dehnes a branched cover of the {x,y) plane, branched 
over the two centres. The inverse image of the x axis {y = 0) consists of the lines 
A = 0 and v = ±7r/2 -|- 2 kn,k G Z. The intersections of these lines, (A,p) = 
(0, ±7r/2 -|- 2k'K) get mapped in an alternating way to the two centres. The inverse 
image of window 3 (the segment between the two centers) corresponds to the line 
A = 0. The inverse image of window 1 corresponds to the lines v = —7r/2 -|- 2 k'K. 
The inverse image of window 2 corresponds to the lines u = tt/ 2+ 2k'K. It is central 
to our analysis that the symbols 1, 2 and 3 are dehned by coordinate lines in the 
separating variables. 

The energy has a single critical point lying on window 3 corresponding to an 
unstable equilibrium lying on the x-axis at the point where the forces exerted by 
the two centres balance. The associated critical value of energy is 

(6) h* =-(^ml-F \/m^)V2 

and serves as a bifurcation value. We view v as an angular coordinate so that 
(A, p) G M X = M/27rZ. Then the Hill region for energy h is the domain 

Hill(h) = {(A, v) : dpx.Pu such that H{X, v^Px^Pu] h) = 0} C M x . 

For h < h* the Hill region consists of two disjoint disks, one disc about each centre. 
These discs merge at the equilibrium for h = h*. For 0 > h > h* the Hill region is 
a connected domain, topologically an annulus x /, wrapping once around the 
cylinder. 

See Fig. 9 and 10 for phase portraits of the two one-degree of freedom Hamil¬ 
tonians for various values of the masses. In these hgures we make a number of 
contour plots in order to indicate how the one-degree-of-freedom systems combine 
to describe the two-degree-of-freedom system. 
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We now collect known facts about the critical values of the one-degree of freedom 
Hamiltonians and the implications these facts have for the critical values of the 
integral map {G,H), see [17] for details. 

Throughout we impose the condition h < 0 since if h > 0 all motions are 
unbounded. We set 


(7) 

and 

( 8 ) 


hx =-{mi + m 2)/2 and =-\mi - m 2 \/ 2 , 

K±a{h) = h + a\mi ± m2|, a = ± 

(mi ± m2)^ 


X±{h) = 


Ah 


Here the first subscript -|- refers to critical values of Hx^ and the — to critical values 
oiH,. 

Critical values of Hx- For h < hx the function Hx{-, ■, h) has exactly one 
critical value, —k++. For hx < h the function has exactly two critical values —k++ 
and -x+. 

Critical values of For h^ < h the function Hi^{-, ■, h) has exactly two 

critical values k and k For h < h^, the function has three critical values k , 

A_+, and x_. 

Critical and Regular values of the integral map. We can deduce 
the critical and regular values of the integral map {G,H) : T*(M x S^) —>■ 
immediately from the above description of the critical values of the one-degree of 
freedom Hamiltionians and the separation of variables. The critical values of the 
map consists of the union of either four or hve analytic curves depending on the 
values of the masses. See figure 11. Two or three of these curves are straight 
lines with slope 1. The remaining two curves are arcs of hyperbolas. The lines 

are g = A++(h) , g = K_+(h), and g = k _(h). The arcs of hyperbolas are given 

by S' = X+{h) for hx < h and g = X-{h) for h* < h < h^. The image of the 
integral map is the region bounded between the leftmost and rightmost of these 
curves (remember: h < 0 always!). This image is divided by the curves into three 
or four simply connected curvilinear domains whose interiors consist of the map’s 
regular values. These regions, henceforth called the “regular regions” are denoted 
by the symbols S’, S, L, P as indicated in the figures. The preimage of a point in 
region S or P is two disjoint 2-dimensional tori, while if the point is in S’ or L this 
preimage is a single torus. The terminology goes back to Charlier [4] and Pauli 
[15]. S and S’ stands for satellite, L for Lemniscate, and P for planetary. 

Our one-degree of freedom systems have period functions Tx{g, h) and T^{g, h) 
defined by an integral over the curve Hx{-, ■, h) = —g and Hu{-, - ^h) = +g. These 
period functions are analytic within each regular region of the {g., h) plane. Our 
tori T{g, h) admit a natural homology basis, or coordinate system, corresponding 
to our separation of variables. We always define the rotation number W in terms 
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Figure 11. Critical values of the energy momentum map {G,H) 
in the symmetric case mi = m 2 = 1/2 (left) and in the asymmetric 
case mi = 3/10, m 2 = 7/10 (right). The regions of regular values 
are marked by S, L, P, for satellite, lemniscate, and planetary type, 
respectively. Blue lines are critical values of H^,] red lines are critical 
values of Hx- 


of this basis. Then our rotation number is 

(9) W{g,h)=mg,h)/Tx{g,h) 

which is a piecewise analytic function. We give explict formulae for W further on. 
One surprise is that W only depends on the value {g, h) of the integral map in the 
case S where there are two tori for a given g, h: the value of the rotation number 
on these two tori is the same. Another surprise is that W is analytic across the 
curve separating the region S from S’. 

See Fig. 12 (bottom) for a plot of W for energy h = —1/4 in the symmetric 
case. For this energy all three types S, L, and P occur. 

The function Tx has discontinuity only along the red singular curves of hgure 11 
so has two ‘branches’ denoted T^s and Txo the hrst having domain formed by the 
union of S, S’ and L, while the second branch has domain P. (In case S’ is empty 
the hrst domain is just S union L.) The function has its only discontinuity along 
the blue singular curves of hgure 11 so also has two branches, denoted T^yo and Tj,r, 
the hrst branch having domain the union of S and S’ and the second branch Ti^r 
having domain the union of P and L. See hgure 13. 

Consequently, the three branches of Ws,l,p of the rotation number are given by: 


Ws{g,h) 


Tuo{g,h) 

Tx3{g,hy 


WL{g,h) 


Tur{g,h) 

Tx3{g,hy 


Wp{g,h) 


Tur{g,h) 

Txo{g,h)' 


For us, the crucial property of the rotation number W is its range. If we hx 
h, and let g vary, the rotation numbers sweep out an interval [Wmin{h), Wmax{h)] 
which have the following dependence on energy h. 
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-0.2 0.0 0.2 0.4 0.6 0.8 1.0 


Figure 12. The period T^, (diverging at 0) and Tx (diverging for 
positive g) (both top) and the rotation number W = T„/Tx (bottom) 
all for h = —1/4, mi = m 2 = 1/2, d = 1. The type of corresponding 
tori is S, L, P from left to right, separated by singularities. 


For type L: 


( 10 ) 


For type S: 


( 11 ) 


h 

min 

max 

Q > h > hx 

0 

00 

hx> h> h* 

Wi(«:++,h) 

00 


h 

^^min 

^^max 

0 > h> h* 

Ws{K_.{h),h) 

00 

h* > h 

1 

1 

Ws{K++{h),h) 


where the critical values fi;++(/i), k_+(/i), k _(h) and h^, hx were defined in equa¬ 

tions (8) and (7) above. 

We now proceed to determine the functional forms of the period functions, and 
hence of the rotation number. We have summarized their forms in Fig. 13. We 
use the standard method of period computation from the classical mechanics of 
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Figure 13. Contours of constant period (top) and Tj, (bottom) 
for the symmetric case mi = m 2 = 1/2 (left) and an asymmetric 
case mi = 3/10, m 2 = 7/10 (right). 


one-degree of freedom systems (e.g. [9]). Let \E'(g,p) denote either Hi^{q,p-,h) or 
Hx{q,p; h). has the form -|- V{q, h). The hrst of Hamilton’s two equations 
reads dq/dr = p so that dr = ^. Solutions must lie on energy curves {T = c}, 
c = ±g, which consist for us of one or two topological circles. Choose a component 
C for a given constant c. The time T to traverse C is the corresponding period T 
which we want to compute. This time is given by the integral 



(One can solve for p in terms of g : p = ±y^2(c — V (g; h) and use the time-reversal 
symmetry to rewrite this as T = 2 =.) In our cases the variable g 

^2{c-V{q;h) ’ ^ 

is either z/ or A and the variable p is either p^ or px The substitution ^ = sin v or 
^ = cosh A converts the integrand y to the integrand ^ where the integral in the 
new variables is around the loop corresponding to our choice C which lies on the 
Riemann surface: 


p2 


2(1 - z^){+g + {mi - m2)z + hz^), 
2(1 - z'^){-g + {mi + m 2 )z + hz'^), 


V case 
A case 


The integral f p over the closed loop is a complete elliptic integral which can be 
expressed in terms of Legendre’s complete elliptic integral K{k^) with modulus k. 
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For the symmetric case mi = m 2 we find 

4 


Tuoig) = 
Turig) = 


4 


ky/—2h 


K(fc^), for g < 0, 
K{l/k^), ioig > 0 , 


where k"^ = 1- y, 

h 

Expressions of elliptic integrals in terms of Legendre’s standard integrals can, e.g., 
be found in [3]. Legendre’s K is a smooth monotonically increasing function that 
maps k"^ G (—cxo, 1 ) to ( 0 , 00 ). 

To write down the period Tx and to treat the asymmetric case mi 7 ^ m 2 intro¬ 
duce 


kl = - + 


g + h 


2 y/ 4 :gh -|- (mi ± 7712)^' 


f±a = 


V 2 


a 


(- 1 - 1 ) 


x/ 4g/i+(mi±m2)2 
h 


where a = 0 or a = 1. Then the A-periods are 

4 


Txsig^h) = 


/+oK(/c+), for g < K++{h) 


Txo{g,h) = 


Vm' 

f+iK{l/kl), for K++{h) <g < x+{h) 


VW\ 

When h > h\ the modulus k\ is monotone in both cases, increasing and decreas¬ 
ing, respectively. Now f+ii{g) is always monotonically increasing, so that Tas is 
monotonically increasing in this energy range. 

Note that since without loss of generality we can set = 1, the Tx periods 

do not change with parameters, only their domain of definition in the integral 
image changes. In particular the periods Tx in the symmetric case mi = m 2 and 
the asymmetric case mi 7 ^ m 2 are given by the same functions. 

For oscillations of u (when mi 7 ^ m 2 ) the period is 

where for h > hi, the domain is 5 ^ < k _|_(h) and the modulus is monotone, while 

for h < hy the domain is g < X- and k'^ passes through zero at h ± (mi — m 2 )- 
The formula is valid in region S and in region S’, in region S it gives the period 
for both tori of type S, around either center. 

For rotations of u (when mi 7 ^ m 2 ) there are two cases depending on whether 
the elliptic curve has 2 real or 4 real branch points. A formula that is valid in both 
cases is defined via 

k^ = - - ^ ^ 

2 2x/{g - h)2 - (mi 


m 2 )^ 
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and 


T^r{g,h) = 


V^{{g — hy — (mi — 


K{k 


When there are fonr real branch points the sign of k‘1 is negative, otherwise k^ 
is positive. Moreover, when h < hi, the modnlus k"^ is monotonically decreasing. 
When h > hy the modulus is monotonically increasing from — oo to 0, which occurs 
at the boundary of the Lemniscate region for g = «:++. 

The functions Tyo{g,h) and Tyj.{g,h) are monotone functions of g, for h hxed, 
because K{k‘^) is monotonically increasing, and k‘^{g) is linear in g and increasing 
for h < 0. So Tyo{g) is monotonically increasing and diverges for g ^ 0-. Now 
l/k'^{g) is monotonically decreasing for > 0 and so as a function of g, Tyr is 
the product of two monotonically decreasing functions, and hence monotonically 
decreasing. It diverges for —)■ 0+. 

Conjecture. The period function T^,^,{g,h) and rotation function W^:{g,h) are 
monotone functions of g within their domains of analyticity. 

Numerical experiments support this conjecture. 

Note that in the S-region the modulus is negative, and approaches —oo at the 
boundary to L. Similarly, in the upper parts of the L- and P-region the modulus 
k"^ becomes negative, and on the boundary between L and S’ it approaches —oo. 

The form of the period functions yield the following facts about the rotation 
number W. For h\ < h < hy the rotation number Wiig) is a monotonically 
decreasing function. At the left boundary of the L region Wl diverges to +oo, 
while at the right boundary of the L region for h > hx the rotation number Wl 
approaches 0. For mi = m 2 the limiting values of W at the other families of 
critical values are as follows, see hgure 15 for graphs of these functions. 

To obtain tables 10 and 11 we need the minimum and maximum of W{-,h) 
within its various intervals of analyticity. For simplicity, in the following formulas 
we have set mi = m 2 = 1/2). We compute that 

• the minimum of Ws{-, h) is ; 


Ws{t^-{h),h) 


+ 1 

2^k(i + ^) 


• the maximum of Ws{-, h) is, if hnite, for h < h*\ 


• the minimum of Wl{-, h) for h* < h < h\ occurs at 
WL{K++{h),h) = -\/-Ah-2 K{-h) 

71 
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Figure 14. Contours of constant rotation number W for the sym¬ 
metric case mi = m 2 = 1/2 (top) and an asymmetric case mi = 
3/10, m 2 = 7/10 (bottom). The spacing of contours is even in W or 
1 /lF, whichever is smaller. 


• the maximum of hFp(-, h) for h\ < h < 0 is: 

In arriving at these facts we used that for the Lemniscate case h) is increas¬ 

ing with g and that Txs^g, h) is decreasing with g, so that Wii^g, h) is decreasing 
with g. We also used that the period diverges at the left boundary of the L- 
region, while the period T\ diverges at the right boundary of the L-region when 
h > hx, hence we arrive at the limiting behaviour oiWi- For a periodic orbit of 
type L or S with rational W = p/q with co-prime p and q the period of the periodic 
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W 



Figure 15. The boundary values of the rotation number W along 
the leftmost curve of critical values (blue, see Fig. 11) which gives the 
minimum in S, and along the rightmost curve of critical values (red, 
see Fig. 11). The three smooth parts of the red curve correspond to 
the rightmost boundary values of tori of type S, L, and P, in that 
order. The red curves give the maximum in S, minimum in L, and 
maximum in P, respectively. Parameters are mi = m 2 = 1/2, d = 1. 

orbit in original coordinates is {pTx + qT^)/2 where the factor 1/2 accounts for the 

double cover. 


8. Finishing up the proof of theorem 2. 

We have established the range of the rotation numbers in the previous section. 
We established almost all of the proof in section 5. All that remains to prove is 
that the windows are situated as in hgure 3. In particular, we have not yet shown 
that the line pairs are spaced 1/2 a unit apart. 

We saw that the two-center system separates into two one degree of freedom 
systems, that the separate torus coordinates correspond to these separate systems, 
with A, Pa corresponding to the the 62 motion and p, pi, to 0i-motion. The windows 
on a given torus are vertical or horizonatl lines in the angle coordinaes since the 
collinear line t = 0 corresponds to A = 0 or p = const and since the axes of 
the torus correspond to one or the other of the separating coordinates: 61 = 
6i{y,py-, g,h), 6*2 = 92{X,p\;g,h). It follows that windows have the form 6*1 = const 
or 62 = const. What constants? In other words, upon being projected onto either 
one-degree of freedom system, the windows correspond to points (or the empty 
set). But which points and how are they separated? 

The explicit maps 61 , 62 as functions of mapping from the regularised coordinates 
{n, A, pu, Pa) are complicated expressions in terms of elliptic functions. Fortunately, 
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this detailed information is not needed to answer our question due to the presence 
discrete symmetries for the separated Hamiltonians. 

The symplectic involution Rx{X,px) = (—A, —px) leaves Hx invariant and hence 
commutes with the time evolution for that flow. The symplectic involution py) 
(vr —z/, —pv) leaves Hy invariant and so commutes with its time evolution (If Rx, Ry 
denotes the extension of these involutions to the full regularized plane: for exam¬ 
ple, Ry{\, h',px,Pu) = (A, n—i',px, —Pv) then the composition of the two involutions 
R = Ry o Rx is an involution on the double cover that leaves the original point 
{x,y,px,Py) hxed.) Rewritten in terms of angle coordinate, 9 = 9i or 92, either 
R must commute with translation (since it commutes with its Hamiltonian flow) 
and thus is of the form R{9) = 9 + p for some constant p. But since each R is a 
nontrivial involution we must have that p = 1 /2. 

For the A motion these window points are 

S’,S,L; the two points (0 ,pa), (0, —Pa) for the symbol 3. exchanged under Rx- 
(Note that Pa = 0 above corresponds to the purely linear motion, and a 
degenerate torus, which we have exculded from considerations.) 

P : No symbol occurs. 

As we just saw, Rx has the form Rx{92) = 02-1-1/2, proving that the corresponding 
window points a half a unit apart, and so the corresponding lines are half a unit 
apart. In other words, the window for the symbol ‘3’ consists of two horizontal 
lines separated by half a lattice unit. 

For the v motion the window points are 

SI: two points of the form (7r/2, ±Pi.) for the symbol 1, exchanged by Ry. 

S2: two points of the form (—7r/2, 9sipy) for the symbol 2, exchanged by Ry. 

P: a point {n/2,py) for symbol 1 and a point {—'K/2,py) for symbol 2, where 
Pj, > 0 and py > 0. There are two more such points with negative momenta 
which belong to a disjoint orbit of Hy (“prograde” or “retrograde” motion) 
The two groups of points are mapped into each other by Ry. 

By an argument identical to the A-case, the two window points in each case are sep¬ 
arated by a half, and so the corresponding vertical lines on the torus are separated 
by half a unit. 

The resulting windows on the flattened torus for the three case L, S, P are thus 
as shown in figure 3. 

The location of the windows relative to the origin of the flattened torus depends 
on the choice of origin of the two angles 9i and 02 and is irrelevant. The inter¬ 
section of windows correspond to collisions, each of which appears twice since the 
regularisation led to a double cover of conhguration space. 

9. Additional Remarks 

Remark: When W = p/q with q even the symbol sequence of type L repeats 
after p + q symbols. When q is odd the 2nd half of the symbol sequence is like the 
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Figure 16. Orbit of type L with IF = 1 and syzygy sequence 1323 
for various values of h = —0.15, —0.2, —0.3, —0.5 from left to right 
and values of mi/m 2 = 0.5, 0.4, 0.3,0.2, 0.1 from top to bottom all 
with mi + m 2 = d = 1. For type L and h > hx = —1/2 all orbits 
exist in such 2-parameter families. 


hrst with symbols 1 and 2 exchanged. For periodic orbits of type S the two halves 
of the symbol sequence are always identical. 

Remark; There are no break-orbits of type L, since for these is always non¬ 
zero. For type S there are break-break orbits for W = p/q with p or q even. When 
both p and q are odd the corresponding symmetric periodic orbit has reflection 
symmetry about the x-axis (but no break points), and the corresponding collision 
orbit is a break-collision orbit. 

Remark; Discrete symmetry of periodic orbits. Let IF = p/q. 

Type L, arbitrary masses; If q is odd, there is an orbit that has a crossing of the 
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x-axis with a right angle. If p is odd, there is an orbit that crosses the origin. In 
addition, if the masses are eqnal and p or q are even, there is an orbit that has a 
crossing of the y-axis with a right angle. 

Type S, arbitrary masses: If p or g are even, there is a break-orbit. If q is odd 
there is an orbit that has a crossing of the x-axis with a right angle. 

Remark: A non-obvious conseqnence of the proof and the description from 
section 7 is that there are no bifnrcations in the orbit strnctnre and the symbol 
seqnences of the L-region when changing the mass ratios or the energy within the 
range {hx, 0). In particniar all periodic orbits persists withont bifnrcations in this 

2- parameter family. This strong resilience of orbits for sufficiently large energy is 
qnite astonishing. It is illustrated for the hh = 1/1 orbit shown for various energies 
and various masses in figure 16. For smaller energies this is not true and orbits 
bifurcate away on the right boundary of region L. 

Remark: Decreasing h from hx to h* (of the saddle equilibrium point) removes 
more and more rotation numbers by increasing the minimum value of the range 
W G (IT™"’, cxd), thus only allowing orbits with more and more consecutive sym¬ 
bols 3, which describe close approaches to the isolated hyperbolic periodic orbit 
above the saddle. It is interesting to note that for any h > h* (even arbitrarily 
close), the whole complexity of orbits generated by {W} passing through [0,1) is 
still present, in fact there are always inhnitely many such intervals. 

Remark: For tori of type P the standard syzygy sequences do not separate 
tori, because symbol 3 does not occur and symbols 1 and 2 simply alternate for 
all allowed rotation numbers W. Instead of insisting on the usual symbols of the 

3- body problem we can turn the construction around and choose a window on the 
torus so that we get similar symbol sequences with a different physical interpreta¬ 
tion. As we already pointed out for separable integrable system the windows are 
best chosen as coordinate lines of the separating coordinate system. We choose an 
ellipsoidal coordinate line A = X*{h) = const such that all type P tori for the given 
energy intersect this line. When approaching the left boundary of the P-region 
the torus shrinks to an ellipse with coshA*(h) = —{mi + m2)/{2h) = hx/h > 1. 
Considering the phase portrait of Hx it is clear that all type P tori intersect the 
corresponding window. However, the intersection points are not spaced half-way 
around the torus, and hence would give a more complicated description. An al¬ 
ternative is to use Pa = 0 which does cut the torus twice with equal spacing. The 
physical meaning of this is to assign a symbol to the tangency with the caustic. 


10. Conclusion 

We have described all syzygy sequences which are realized in the two centre 
problem. They are encoded by Sturmian words as per theorem 1 and 2. 
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11. Some Open Problems. 

Orbit Counts. This work began as a warm-up for understanding the syzygy 
sequences arising in the full planar three-body problems. Does it shed light on 
that problem, or on restricted versions of it? We take a peek at this question from 
the perspective of orbit counting. 

Let N{L) be the number of periodic orbits of length less than L in a dynamical 
system. For chaotic systems such as Axiom A systems N{L) grows exponentially 
with L, while for integrable systems one expects N{L) to grow polynomially. Var¬ 
ious theorems have been established relating N{L) to the topological entropy of 
a system. In planar three-body problems we will instead let N{L) be the num¬ 
ber of distinct syzygy sequences of length L which are suffered by some periodic 
orbit during one period. (The original N{L), the number of orbits themselves 
is uncountably inhnite since the orbits of a particular topological type, forms a 
continuum in an integrable system, sweeping out all of a torus, or a torus minus 
collisions. So we do not want to count individual orbits.) Being integrable, we 
expect that N{L) grows polynomially for the two centre problem. Indeed, N{L) 
grows like L^. To estimate this growth, recall that altogether P orbits yield a 
single syzygy sequence 1212 .... To count the number of sequences due to L and S 
orbits, recall that these sequences are encoded by their rational rotation number 
W = p/q, W > 0 with the length L of the sequence being 2(p -|- q). So N{L) is 
bounded by the number of lattice points (p, q) inside the diamond |a;| -1- \y\ < 
which is We count both the S and L type, yielding the bound N{L) < ^L‘^ + 1. 

On the other hand, if we use the strong force potential —mi/r\ — m 2 /r\, then 
one can prove using variationally methods, that all syzygy sequences are realized 
in the two centre problem except for those with “stutters”, meaning two of the 
same type of symbols adjacent to each other. One estimates N{L) > 3 * 2^“^ for 
the number of such symbols. 

The two-centre problem hts within various one-parameter families of 2-degree 
of freedom Hamiltonian systems Hx(x,y,Px,Py)- For example, by “spinning” the 
primaries we can £t it into a one-parameter family, parameterized by the spin 
rate of the primaries, which at the other endpoint becomes the restricted three- 
body problem. Or, by changing the potential, we can £t the problem into a 
one-parameter family which includes strong-force two center problem. How does 
Nx{L) depend on A for these other problems? 

Sturn and twist As remarked above, Sturmian words arise generically in 
twist maps [2]. This “Sturmian thread” connecting the two centre problem to 
general twist maps gives us hope that enough twist-map structure persists as we 
turn on the spin parameter so that many of the periodic orbits or quasi-periodic 
irrational Sturmian type orbits will persist all the way to the restricted three body 
problem. 
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